Finite-size and Particle-number Effects in an Ultracold Fermi Gas at Unitarity 
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We investigate an ultracold Fermi gas at unitarity confined in a periodic box V = L 3 using 
renormalization group (RG) techniques. Within this approach we can quantitatively assess the long 
range bosonic order parameter fluctuations which dominate finite-size effects. We determine the 
finite-size and particle-number dependence of universal quantities, such as the Bertsch parameter and 
the fermion gap. Moreover, we analyze how these universal observables respond to the variation of 
an external pairing source. Our results indicate that the Bertsch parameter saturates rather quickly 
to its value in the thermodynamic limit as a function of increasing box size. On the other hand, we 
observe that the fermion gap shows a significantly stronger dependence on the box size, in particular 
for small values of the pairing source. Our results may contribute to a better understanding of finite- 
size and particle-number effects present in Monte-Carlo simulations of ultracold Fermi gases. 



I. INTRODUCTION 

Ultracold atoms provide accessible and controllable 
systems to study quantum many-body effects pQ. Of- 
ten a single parameter, the s-wave scattering length a s , 
describes the microscopic interactions completely, and 
can be tuned by means of an external magnetic field in 
the presence of a Feshbach resonance. In the context 
of ultracold fermions [2], this gives rise to the BCS-BEC 
crossover [3] , which interpolates between the cornerstones 
of quantum condensation phenomena, namely fermion 
superfluidity via Cooper pairing and Bosc-Einstein con- 
densation of molecular bosons, and has been probed in 
milestone experiments [3]. 

The intermediate regime of the BEC-BCS crossover, 
where the scattering length is large and no obvious small 
expansion parameter exists, represents a major challenge 
for theoretical many-body approaches. Recently, experi- 
ments in this so-called unitary regime are reaching a level 
of precision [5] that calls for an improved quantitative 
understanding on the theory side. 

Quantitative theoretical understanding can be ob- 
tained from numerical (Quantum) Monte Carlo simula- 
tions with the advantage of not relying on specific ap- 
proximation schemes, see e. g. Refs. [6-12J. On the other 
hand, various functional techniques exist which allow to 
associate specific physical mechanisms with the result- 
ing numbers for universal observables. With these tech- 
niques, the complete phase diagram has been studied. 
In particular, e-expansions [HflfTB"] . 1/iV-expansions [TTJ 
118] . i-matrix approaches |19j . Dyson-Schwinger equa- 
tions HOI HI], 2-particle irreducible methods [52], and 
renormalization-group flow equations [2"5r[5T] have been 
employed. 

Despite these tremendous theoretical efforts, our un- 
derstanding of the limit of large scattering length is still 
not complete. In the following we shall focus on this 
limit for a spin-balanced two-component Fermi gas at 
zero temperature, and study the effects of both a finite 
volume and a finite particle number on universal quanti- 



ties, such as the Bertsch parameter and the fermion gap. 
A key motivation for such an analysis is to foster the 
comparability of analytical approaches to data from lat- 
tice (Monte-Carlo) simulations which are performed in a 
finite volume. In the context of resonantly interacting 
Fermi gases, finite-size and particle-number effects are 
currently under investigation using different Monte-Carlo 
approaches [TTJ [T2] I52"H55] . Of course, a comprehension 
of such effects may also provide useful information for 
an analysis of data from Quantum Monte-Carlo stud- 
ies [3B] of the Tan relations [37]. Let us also mention 
that finite-volume effects in atomic few-body systems 
have been studied recently using an effective field the- 
ory approach [38Tl40j . Finally, we would like to add that 
a study of finite-size and particle-number effects arising 
from a given experimental setup may eventually allow us 
to make better contact between theory and experiment, 
see e. g. Refs. [HI [42]. 

For our analysis we employ a functional renormaliza- 
tion group (FRG) approach and apply techniques similar 
to those used for the analysis of finite-volume effects in 
QCD low-energy models [4"5H4"T] . For our purposes, such 
an approach is advantageous since it allows us to inves- 
tigate the "transition" between the finite system and its 
continuum limit. 

Within our RG approach, it is possible to continuously 
follow the change of relevant degrees of freedom by grad- 
ually integrating out the high momentum modes. Here, 
we are interested in quantifying long wavelength physics 
on the scale set by the finite size L, which typically ex- 
ceeds the length scale set by the inverse Fermi momen- 
tum kp , or by the chemical potential /i defining the Fermi 
surface for positive values, yfjl < fcp, which is bounded 
from above by the Fermi momentum. The physics on 
these momentum scales well below the Fermi surface is 
dominated by bosonic fluctuations, which at zero temper- 
atures physically are fluctuations of the superfluid order 
parameter in the crossover problem. 

We shall see below that our analysis fully supports 
this basic picture: Indeed, we find that the effects of 
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a finite system size are rather weak in a mean field ap- 
proximation, which takes only fermion fluctuations into 
account but omits those of the order parameter. Tak- 
ing the bosonic fluctuations into account, we find that 
the finite size effects on the Bertsch parameter are still 
rather weak, whereas the fermion gap is more sensitive. 
This complies with the general result that long wave- 
length phase fluctuations may have a strong impact on 
the condensate. In particular, as is well known, in lower 
spatial dimensions, these fluctuations may indeed com- 
pletely destroy the long range order. 

From this argumentation, we see that our investiga- 
tion can indeed help to quantitatively estimate the finite 
size errors of Quantum Monte Carlo simulations, which is 
notoriously difficult to improve on in lattice simulations. 
In particular, both the precise scaling of an observable 
and the value of its relative deviation for a given fixed 
size from the continuum results may be valuable in view 
of direct comparison, even if the absolute value of this 
observable in the continuum limit may not be captured 
with high precision by our analytical approach due to an 
incomplete treatment of the short-range physics. 

As an important technical development in our study, 
we allow for a finite external (pairing) source J which 
couples to the order-parameter field. The relevance of 
such an external source for studies in a finite volume has 
been pointed out in Refs. [32l [48] . In any case, the in- 
clusion of a pairing source allows us to control symmetry 
breaking in a finite volume. As done in lattice (Monte- 
Carlo) simulations, see e. g. Ref. [3], we shall choose the 
boundary conditions of the fermions in the spatial direc- 
tions to be periodic. Moreover, we consider the grand 
canonical ensemble for our studies. This implies that our 
choice for the chemical potential fixes the (average) par- 
ticle number of the system. 

Finally, we add that we do not study possible finite- 
density corrections in lattice simulations which arise from 
the discretization of space-time and are associated with 
modifications of the (continuum) dispersion relation on 
the lattice. This type of finite-size effects has been 
addressed in a recent Dynamical Mean Field (DMFT) 
study [49] . 

The paper is organized as follows: In Sect. [TTJ we 
discuss general aspects of finite-size effects in ultracold 
Fermi gases at unitarity. A mean-field analysis of finite- 
size and particle-number effects is then presented in 
Sect. |III| In Sect. |IV] we study the impact of order- 
parameter fluctuations in a finite system, i. e. we take 
into account effects arising from the Nambu-Goldstone 
modes in the spectrum of the theory. Our conclusions 
and a brief outlook are given in Sect. [V] 



II. GENERAL ASPECTS OF EFFECTS OF A 
FINITE VOLUME AND A FINITE EXTERNAL 
SOURCE 

Let us start our general discussion of finite-size and 



particle-number effects in a Fermi gas at unitarity by 
considering a standard mean-field approximation. We 
emphasize, however, that the general finite-size formal- 
ism developed here is directly applicable to the analysis 
beyond mean field in the Sect. IV The mean- field ap- 



proximation amounts to dropping contributions associ- 
ated with fluctuations of the order-parameter field. In 
Gross-Neveu-type models, for example, such an approxi- 
mation corresponds to the leading order contribution in 
a systematic expansion in powers of 1/Nf, where Nf de- 
notes the number of fermion species |50j . Such an ap- 
proximation already allows us to discuss finite-size effects 
as far as they relate to the fermionic modes of the sys- 
tem. The role of the dynamics of bosonic bound-states 
of fermions will be discussed quantitatively in Sect. |PV| 

The microscopic (classical) action S of a resonantly 
interacting Fermi gas can be written as follows: 

S[^, ip] = J dr J d 3 x {v f (d T - V 2 - /x) i/> 

+ *M^)W>V)} »(!) 

where ip T = {tjj*, ipi) an d A^ denotes the bare four- 
fermion coupling. The chemical potential is given by /z. 
Note that 2m = 1 in our conventions, where m is the 
mass of the fermions. 

The dimensionless renormalizcd four-fcrmion cou- 
pling A,/, ~ \ipk in Eq. ([!]) is related to the s-wave scat- 
tering length a s , where k is the RG scale. To be specific, 
we have 



A,/, = 



8ttA 



-rcg 



A 



(2) 



in the limit of a broad Feshbach resonance, where A de- 
notes the ultraviolet (UV) cutoff. The constant c reg , > 
depends on the employed regularization scheme. For ex- 
ample, we have c rcg . = 2/tt for the sharp UV cutoff. 

In order to study the RG flow of the quantum effective 
action T, we employ a non-perturbative RG equation, 
the Wetterich equation |51j . The effective action then 
depends on the RG scale k (infrared cutoff scale) which 
determines the RG "time" t = ln(fe/A). For reviews on 
and introductions to this functional RG approach we re- 
fer the reader to Refs. [2U1 \52lWZ\ . 

For our analysis of finite-size and particle-number ef- 
fects, we employ the following ansatz for the effective 
action: 

r k [tf,1>,<p] = JdrJ d 3 x{^ (fl r -f 2 -/i)^+ro*pV 



-K [(p*^tpi-ipip^ipl\ } , (3) 



where J denotes an external pairing source and the coef- 
ficient of the cubic fermion-boson interconversion term, 
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hp, is a Feshbach or Yukawa coupling. Note that the 
source term oc J explicitly breaks the underlying U(l) 
symmetry of the microscopic theory. Formally, this 
ansatz is inspired by the partially bosonized formulation 
of the action ([!]) , which can be obtained by introducing 
a complex scalar field ip ~ (h ip /fh 2 )(tpt' t Pl) m t° the path 
integral by means of a Hubbard-Stratonovich transfor- 
mation. Physically, this model describes the physics of 
fermionic atoms close to a Feshbach resonance in an in- 
tuitive way, where the to^, plays the role of the detuning 
from a bosonic bound state level, described by the com- 
plex field ip, at the microscopic scale. The bare couplings 
h v and fh v are a priori constants at our disposal and 
are chosen such that the four-fermion interaction term in 
Eq. ([lj vanishes identically, = —h^L/fh 2 . The model 
becomes physically equivalent to the model ([!]) in the 
broad resonance limit, where h v — > oo, while the ratio 
h^/rh 2 , is kept fixed [20] on the microscopic scale, which 
we consider in this work. For our purposes below, it is 
convenient to consider two real- valued fields ip± and ipi 
instead of the two complex- valued fields ip and ip* , where 
(p = (ipi+itp 2 )/ \/2. We would like to add that a partially 
bosonized formulation is beneficial since it allows us to re- 
solve (parts of) the momentum dependence of fermionic 
interactions by means of a straightforwardly accessible 
derivative expansion of the effective action. 

Due to quantum fluctuations, the cubic interaction 
generates kinetic terms for the bosonic fields in the RG 
flow, 

Zyd T ip and Z>*VV, (4) 

even if these terms have been set to zero at the initial 
RG scale A. For example, taking these fluctuations into 
account by following the RG flow down to k -> 0, and 
concentrating on the BEC regime of the crossover for a 
moment, where the bosonic field describes tightly bound 
molecular states, their ratio is given by Z^jz\ — > 1/2 or, 
with dimensions restored, l/(4m). Indeed, this signals 
propagating bosonic bound states of double mass of the 
fermionic atoms. More generally, the leading momentum- 
dependence of the four-fermion vertex in Eq. ([I]) in the 
pairing channel is encoded in these kinetic terms. In lead- 
ing order in a systematic expansion of T in derivatives, 
the RG flows of the wave-function renormalizations zj^' 
are essentially determined by a purely fermionic one- 
particle irreducible (1PI) diagram in the U(l) symmetric 
regime. Thus, we have 

d t ]nZl< ± = -c v h 2 v , (5) 

where c v is a positive constant |29) . Contributions from 
1PI diagrams with internal boson lines are suppressed 
in the symmetric regime due to the large (renormalized) 
boson mass parameter m v = m^/(Z^k 2 ), at least in the 
limit of a broad Feshbach resonance. For a similar reason, 
the running of the fermionic wave-function renormaliza- 
tions is subleading since corrections due to 1PI diagrams 
with at least one internal boson and fermion line are sup- 
pressed, both in the U(l) symmetric regime as well as in 



the regime with broken U(l) symmetry in the ground 
state. Finally, the RG flow equation of the Yukawa cou- 
pling also assumes a simple form in the symmetric regime, 
since it is only driven by the anomalous dimensions of the 
bosonic fields: 

d t h% = (r£ - l)h% , (6) 

where h 2 = h%,/(Z£k), see Ref. [2H]. Note that we are 

free to choose either z\ or to renormalize the Yukawa 
coupling. In the present work, we only take into account 
the running of z\ and set z\ = Z^. We would like to 
add that the values of low-energy observables, such as 
the fermion gap and the Bertsch parameter, do not de- 
pend on the RG flow of the wave-function renormaliza- 
tions in the mean-field limit. Beyond this limit, however, 
the running of the latter affect the running of physical 
observables. 

The universal behavior associated with the limit of a 
large (infinite) s-wave scattering length a s is linked to 
the existence of a non-trivial UV attractive fixed-point of 
the theory [TTI 125] . In a given regularization scheme, the 
value of the four-fermion coupling given in Eq. (§ can be 
associated with the ratio of the UV fixed-point values h* v 
and to* of the Yukawa coupling and the mass parameter, 
respectively. Recall that — —h^/rh 2 ,. In order to 
study the limit of infinite s-wave scattering length, we 
therefore have to choose the initial conditions for the RG 
flow equations of rh v and h v such that they are close to 
their UV fixed point values, see e. g. Refs. [2"51 l2"9l R)2"] . 

In our ansatz ^ we also allow for a term cx A^ which 
describes four-boson interactions. At the initial RG 
scale A we set the coupling \ v to zero. For scales k < A, 
this coupling is then generated due to quantum correc- 
tions, even if we do not take into account fluctuations of 
the bosonic fields. Boson self-interactions of higher or- 
der (ip*ip) m with m > 2 are also generated. We add 
that the bosonic self-interaction terms (with m > 1) 
essentially parameterize the effective potential (order- 
parameter potential) U. In the following we drop the 
contributions from higher-order interactions (to > 2) 
for simplicity. We rush to add that the RG flows of 
the associated couplings are decoupled in the mean-field 
limit in the U(l) symmetric regime at scales fc > fcsB, 
where ksB denotes the scale associated with spontaneous 
U(l) symmetry breaking. 1 Thus, the determination of 
the scale fcgs does not suffer from this approximation, 
at least in the mean-field limit. However, the (exact) 
value of a given low-energy observable O ~ fcg B may re- 
ceive contributions from boson self-interactions of higher 



For a comprehensive discussion we refer the reader to Ref. 1501 
where this issue has been studied in detail in the Gross-Neveu 
model. Since these arguments only rely on very general proper- 
ties of partially bosonized formulations of fermionic theories, the 
same reasoning also holds in the present study. 
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order. 2 As we shall see below, these corrections to the 
fermion gap and the Bertsch parameter are small. Again, 
this holds at least in the mean-field approximation. 

As already indicated above, we also include a source J 
for the field ip± in our RG study. This term breaks explic- 
itly the underlying U(l) symmetry of the theory. How- 
ever, such a linear symmetry breaking term remains un- 
changed in the RG flow [53]. Instead of including such a 
term in the RG flow, one may therefore simply study the 
evolution of the effective action without a finite external 
source J. Explicit symmetry breaking can then be taken 
into account after the quantum fluctuations have been in- 
tegrated out on all scales. This strategy has been followed 
in the context of low-energy models of QCD, see e. g. 
Rcfs. 64, 65 . In particular, such an approach is perfectly 
suited when one has access to the full effective poten- 
tial U{p) = TT[p]/L 3 including bosonic self-interactions 
of arbitrarily high orders, where p — {(p*(p); Cp denotes 
a spatially constant background field and 1/T is the ex- 
tent of the system in Euclidean time direction. Using 
a low-order expansion of the effective potential as indi- 
cated in Eq. ^ , such an approach still yields reasonable 
results in an infinite volume and for small values of the 
source J. In a finite- volume study, however, the situation 
is different. Without a finite value for J, fluctuations of 
massless Nambu-Goldstone bosons associated with spon- 
taneous U(l) symmetry breaking restore the symmetry 
in the (deep) IR limit. This implies that there is no 
spontaneous U(l) symmetry breaking in a finite volume. 
Therefore we include a finite source J on all scales in the 
RG flow to control U(l) symmetry breaking in a finite 
system. 3 This renders the order-parameter (<f>)j ~ (4>i)j 
finite on all scales. 4 

Let us briefly discuss a subtlety concerning the ex- 
pectation values {(f)) j and {\(f)\) = {\(f>\)j=o. In the 
limit J — > 0, the minima of the order-parameter poten- 
tial U are degenerate and describe a circle of radius {\<f>\) 
in the plane spanned by <p\ and <p2- For a finite value of J, 
the degeneracy of the ground states is completely lifted 
and the mass mc of the pseudo Nambu-Goldstone modes 
is finite. We will now argue that \{(f>)\ := \ \\mj^ {(f)) j\ 
and {\4>\) agree in the infinite-volume limit. In this limit, 
both quantities may therefore be considered as equivalent 
order parameters for spontaneous U(l) symmetry break- 
ing. In order to understand better the relation between 
these two quantities, let us first consider a system in a 
finite periodic volume V. In this case, loosely speaking, 



2 Here, we have tacitly assumed that O has the dimension of en- 
ergy. 

3 Due to our definition of J, we have (<j>) j = (<j>i)j/\/2 > 
and (<j) 2 ) = for J < 0. 

4 In this work, we employ the convention that the field rf> used 
in the expectation value (...) denotes the quantum field which 
appears in the path integral. This quantum field should not be 
confused with the so-called classical field ip which appears in the 
effective action. 



one finds that the fluctuations of the spatial (momentum) 
zero-modes of the Nambu-Goldstone fields along the cir- 
cle of energetically equivalent ground states tend to re- 
store the symmetry. Schematically, these fluctuations are 
suppressed according to 5 ~ exp(-(mG\ / ^) c ), where ttiq 
denotes the mass of the (pseudo) Nambu-Goldstone par- 
ticles for J ^ and c > 1 is a positive constant [66j . 
Considering the limit toq ( J 0) before taking the 
limit V — > oo, we then anticipate that {4>) j-q averages to 
zero since the fluctuations of the zero modes are not sup- 
pressed. On the other hand, these fluctuations of the zero 
modes are projected out in the computation of {\<f>\). In 
a finite volume at sufficiently low temperatures, {\(f>\) can 
therefore be finite for J = 0, even if we have (</>)j = o = 0. 
In fact, we always have {(f)) j—q = in a finite volume. 
Considering now the limit V — ¥ oo and then uiq — > 0, 
fluctuations of the Nambu-Goldstone along the circle of 
energetically equivalent ground states are exponentially 
suppressed and {4>)j=o remains finite. In particular, we 
have \{<f>)\ = (|0D- For sufficiently large volumes, how- 
ever, it is reasonable to assume that the deviations of 
{\(j)\) from the value of \{4>)\ are small in the thermody- 
namic limit. 6 In the following we always use {<f>)j for 
J 7^ as a (pseudo) order parameter for U(l) symmetry 
breaking in a finite volume. Note that the fermion gap 
and the ground-state energy of the Fermi gas depend on 
the actual value of {(f)). j- 

Finally a comment on inhomogeneous ground states 
is in order. In a finite volume and for a finite particle 
number the ground state of the theory might be inhomo- 
geneous, implying that {4>)j has a non-trivial dependence 
on the spatial coordinates. 7 From studies of relativistic 
theories, there is indeed direct evidence that finite-size 
effects may alter the phase structure of a given theory. 
For example, lattice studies of the 1 + Id Gross-Neveu 
model show that the finite-temperature phase diagram 
of the uniform system is modified significantly due to the 
non-commensurability of the spatial lattice size with in- 
trinsic length scale of the inhomogeneous condensate [6 9) . 
However, a study of the emergence of an inhomogeneous 
ground state as a function of N and L is beyond the 
scope of the present work. From here on, we shall as- 
sume that the ground state of the theory is homogeneous, 
even for finite values of N and L. Thus, we consider the 
(pseudo) order-parameter {(f)) j to be independent of the 
space-like coordinates. This amounts to exploring only 



5 For the Goldstone bosons with linear dispersion, the relativistic 
power counting, where mQ corresponds to an inverse length, is 
appropriate. 

6 Note that fluctuations of the Nambu-Goldstone are not included 
in a mean- field study. This implies that (<j>) j remains finite in a 
finite volume at sufficiently low temperatures, even for J — > 0. 

7 The ground state of a resonantly interacting Fermi gas may be 
inhomogeneous even in the infinite-volume limit. For example, a 
Sarma-type 1671 or Fulde-Ferrell-Larkin-Ovchinnikov 68 phase 
may exist for unequal chemical potentials of the spin-up and 
spin-down fcrmions. 
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leading-order effects on physical observables arising from 
the presence of a finite volume V and a finite particle 
number N. At least close to the thermodynamic limit, 
it still seems reasonable to expect that corrections due 
to possible inhomogeneities of the ground state are sub- 
leading. 



III. MEAN-FIELD ANALYSIS OF FINITE-SIZE 
AND PARTICLE-NUMBER EFFECTS 

Let us now discuss our RG study in the mean-field limit 
in more detail. The flow equation of the order-parameter 
potential in the mean-field limit can be derived along 
the lines of Ref. [21] , where the thermodynamic limit has 
been studied in detail. For a system in a finite volume 
we then obtain: 

d t U(p, J,L,/i) = -2k 5 (B> + B<)s F , (7) 

where the functions Bp and Bp are given by 

1 , 



({kL) 2 ~{2v) 2 q 2 + pL 2 ) 

x6((27T) 2 q 2 -^L 2 ) , (8) 



x0 (fiL 2 - (2tt)V 2 ) • (9) 

These functions count the momentum modes. Recall that 
we use periodic boundary conditions for the fermions in 
spatial directions. Furthermore, we have to specify the 
function s F at T — 0, 



fc 4 + Kp 



(10) 



The fermion gap is given by A = h^po — ft. 2 \((j>)j\ 2 . To 
regularize the theory, we use here and in Sect. |IV| an opti- 
mized regulator function |29j . For details on optimization 
of RG flows, we refer the reader to Refs. [5711701172] . 

A quantity of utmost importance in our study is the 
fermion density n. The RG flow of the density n can be 
deduced from the flow equation of the effective potential. 
For details, we refer the reader to Ref. The (average) 
particle number N is then given by N = nL 3 . The initial 
condition for the RG flow of the density is given by the 
one of the free Fermi gas nf reG , which is determined by 
our choice for the chemical potential. To be specific, nf reo 
is given by 



L 3 



8 (fiL 2 -(2ir) 2 q 



2\ i^L'-^ca) 



3tt 2 



(11) 



where (ft G Z and s is the spin-degeneracy factor. For a 
two-component Fermi gas, we have s — 2. Apparently, 



nfree is discontinuous for finite values of L. We shall come 
back to this below. 

For clarity, let us discuss the definition of the Bertsch 
parameter and the (dimensionless) fermion gap which are 
at the heart of our interest in the present work. In the 
limit of large s-wave scattering lengths a s , (dimension- 
less) IR observables are universal, i. e. they do not de- 
pend on the actual value of the density n, provided that 
we consider J — > in the thermodynamic limit. In our 
numerical studies, /j, represents a dimensionful parameter 
which can be adjusted by hand. For a given value of /i, 
we then compute the density n of the interacting system. 
This density can then be related to the Fermi energy e F 
of a free Fermi gas with the same density: e F = kp, where 
the Fermi momentum is given by k F = (Z^n) 1 / 3 . bmce 
we consider the limit a s — > oo, the energy E of the inter- 
acting system is proportional to /i, as it is the case for the 
free system. Thus, we have E/N — £Ep/N, where £ is 
the so-called Bertsch parameter. For the free system, we 
have E F = 4 /2 V7(57r 2 ). The energy per particle is given 
by Ep/N — (3/5)e F . Since \i ^ e F for the interacting 
system, it follows that 



E 
N 



:£e F and fj, = £e F = £(3?r 2 n) 3 , (12) 



where E/N :— (3/5) ji. The dimensionless fermion 
gap A can be defined accordingly by using the Fermi 
energy of the free system as a reference scale, A = A/e F . 
In our studies of a Fermi gas in a finite volume V = L 3 in 
the presence of a finite pairing source J, the energy E 
depends on three scales, namely p, L and J. For a 
given set of values for /i, L and J, we can then com- 
pute the density n = n(fi,L,J) and the average par- 
ticle number N = nL 3 of the interacting system. In 
order to "measure" the energy of the interacting Fermi 
gas, we again use the Fermi energy e F of a free uni- 
form Fermi g reference scale. For convenience, we 
define the energy per particle of the interacting system 
as E/N := (3/5)/x(iV, L, J). In analogy to the thermo- 
dynamic limit, we can then define the Bertsch parameter 
as follows 



E(N,L,J) 3^ Ar T TX 

AT - -£(N,L,J)e F and 
iv 5 

»{N,L,J) = £(N,L, J)(37r 2 n)i , 



(13) 



where n is the density of the interacting finite system. 
From £(N, J, L) we can then read off how the energy 
(density) of the free uniform system changes when we 
vary N, J and L. The dimensionless gap in a finite vol- 
ume is again defined by using the Fermi energy e F = 
(37r 2 n) 2 / 3 of a free uniform g clS ctS Hi reference scale, where 
the density of this free gas is identical to the density 
of the interacting finite system. Our observations re- 
garding the volume dependence of the Bertsch param- 
eter and the fermion gap can in principle be used to 
guide Monte-Carlo simulations of the type discussed in 
Ref. [48] . However, our observations may also provide 
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FIG. 1: Relative shift of the Bertsch parameter 5£ as a 
function of Snj = (n.j — rij=o)/nj=o in the thermodynamic 
limit. The black solid line depicts the numerical results from 
our mean-field study for Jfcp 7//2 = 0.007... 0.7 (or equiva- 
lently J^~ 7/4 = 0.003 ... 0.3). The triangles are the results 
for Snj from our study including order-para met er fluctua- 
tions, where Jkl /2 = 0.009 . . . 1.023, see Sect. 



IV 



While the 



absolute values for £ differ from the ones of the mean-field 
study, the results for <5£ are indistinguishable from the mean- 
field results on the scale of the plot. The red dotted line and 
the blue dashed line represent the results from the expansion 
of 5£ up to first and second order in Snj, respectively, see 
Eq. (fl4|. 



useful insights into the N- and L-dependence of results 
from other Monte-Carlo approaches. 

Before we study the Fermi gas in a finite volume, we 
briefly discuss the effect of a finite (pairing) source J in 
the thermodynamic limit. To this end, we analyze the 
J-dependence of the fermion gap and the Bertsch pa- 
rameter. While the Bertsch parameter can be viewed as 
a thermodynamic observable, which depends only indi- 
rectly on J, the fermion gap plays the role of the order 
parameter for U(l) symmetry breaking and is therefore 
directly related to the pairing source J. For J — > 0, we 
obtain £ rs 0.60 and A/ep rj 0.63 within our mean-field 
approximation. 8 At this point, we would like to remind 
the reader that we are not aiming at an exact determina- 
tion of these quantities in the present work but only at 



Our result for the fermion gap does not agree with the "exact" 
mean-field result for this quantity, see e. g. Refs. [3] 1281 I73| . 
since we have dropped contributions from bosonic self interac- 
tions {<p*<p) m of higher order (m > 2), see our discussion above. 
In the regime with broken U(l) symmetry, the RG flows of 
bosonic self-interactions (ip*ip) m are coupled even in the mean- 
field limit due to the presence of a finite vacuum expectation 
value (cj>). Note that the higher-order interactions (m > 2) are 
included in a standard mean-field approach. In contrast to the 
fermion gap, however, the value of the Bertsch parameter ap- 
pears to be less sensitive to the inclusion of such higher-order 
interactions. 



FIG. 2: Relative shift of the fermion gap SA as a func- 
tion of Snj — (n.j — nj=o)/nj=o in the thermodynamic 
limit. The black solid line depicts the numerical results from 
our mean-field study for Ju 7 ^ 4 — 0.003 . . . 0.3 (or equiva- 
lently Jklp = 0.007 . . . 0.7.) The black dotted line represents 
the result of the expansion of SA up to first order in Snj, 
where the expansion coefficient has been obtained from a fit 
to the numerical data. The corresponding results from an ap- 
proximation beyond the mean-field limit are depicted by the 
cyan-colored (light gray) lines. 



an understanding of the scaling behavior of these quan- 
tities associated with a variation of J and L. For finite 
values of J, we have two scales in our theory, namely k-p 
and J. The pairing source J yields a contribution to 
the density and therefore hp depends implicitly on J. In 
addition, J may contribute directly to the fermion gap, 
i. e. A = A(fcp, J). As a consequence, the dimensionless 
gap A/ep is no longer a constant but depends on the ac- 
tual values of the source J and the Fermi momentum k-p. 



For (small) finite values of J, expanding Eq. ( 12 ) we find 
that the relative shift of Bertsch parameter o£takes the 
form 

5H = = -\{8n,) + l{8n,Y + ((Snj) 3 ) , (14) 

where £ = O=0 and Snj = (nj - nj =0 )/n J=0 . Here, 
nj denotes the density in the presence of the source J. 
Recall that the Bertsch parameter is independent of fcp 
for J — > 0. Thus, the variation of St; can be expressed 
solely in terms of the shift of the density Snj which de- 
pends on both fcp and J. We add that Snj tends to zero 
for increasing fcp, if we keep the source J fixed. In fact, 
an increase of fcp for fixed J corresponds to an increase 
of the total density while effectively keeping the contri- 
bution from the source J fixed. On the other hand, Snj 
becomes larger for increasing J and fixed fcp. 

The expansion coefficients in Eq. ( 14 ) are exact, i.e. the 



general form of the expansion does not rely on the mean- 
field approximation. In fact, the expansion coefficients 
are independent of the approximation scheme. In prin- 
ciple, Snj can be also expanded in powers of J. For an 
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analysis of lattice data, however, an expansion in powers 
of Snj seems to be more suitable. In Fig. [I] we show 6£ as 
a function of Snj as obtained from our mean-field approx- 
imation (black/solid line). The linear approximation and 
the quadratic approximation are given by the red/dotted 
line and the blue/dashed line, respectively. We observe 
that the linear approximation is justified for Snj < 0.05. 
For Snj > 0.05, higher order corrections become impor- 
tant. The knowledge of the range of validity of the linear 
approximation might be useful for the analysis of data 
from lattice simulations. In these simulations only data 
for finite values of J might be available such that the 
limit J — > needs to be extracted from a fit to the avail- 



able data obtained for J ^ 0. The expansion (14) may 



be well-suited for such a fit since only nj =0 and £o enter 
as free (fit) parameters. 

For the dimensionless fermion gap A = A/ep one 
might be tempted to consider an expansion in powers 



of the (dimensionless) source j = Jk } 



-7/2 



Such an ex- 



pansion seems to be natural for this quantity since the 
source J contributes directly to the fermion gap A. This 
can be immediately seen from the effective action ^ and 
the fact that tp± ~ V'tV'i- From a practical point of view, 
however, it seems to be more appropriate to also consider 
an expansion in powers of Snj rather than j: 



SA 



-A 



A 



S^ (Snj) + 0((Sn d 



(15) 



where the S^ denotes a dimensionless expansion coef- 
ficient. In Fig. [2] we present our results for the relative 
shift of the gap SA as a function of Snj. From a lin- 
ear fit to our mean-field data, 9 we find S^ « 0.749, 
see red/dotted line in Fig. [2j A (linear) fit to our results 
from a study beyond the mean-field limit 10 (see Sect. IV) 



yields S^ « 1.563. We conclude that the expansion co- 
efficients depend on the truncation scheme. In fact, our 
results suggest that long-range order-parameter fluctua- 
tions absent in a mean-field study tend to increase the 
slope 8^ . An exact determination of this coefficient is 
beyond the scope of this work. We only state that the 
linear approximation is in good agreement with the nu- 
merical data for Snj < 0.4 in both cases. 

Now that we have clarified the role of the source J in 
the thermodynamic limit, we turn to a discussion of the 
effects of a finite volume (in the presence of the source J) . 
Since we consider a Fermi gas in a volume V = L , 
the momenta of the fermions are discrete. These dis- 
crete momenta define a lattice of equidistant sites in 
momentum space. A specific momentum then corre- 
sponds to a specific site of this lattice. In Fig. [3] we 
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FIG. 3: Number of spin-up and spin-down fermions N 
2N t = 2iV; of a free Fermi function of fj, L 2 / (2tv) 2 . 
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of a non-interacting gas as a function of the dimension- 
less quantity uL 2 /(2ir) 2 . Increasing \xL 2 , we find that a 
jump/discontinuity occurs in the particle number N. At 
such a discontinuity the number of modes contributing to 
the partition function increase. In other words, the num- 
ber of modes enclosed by the Fermi sphere with radius k-p 
changes only when \xL 2 assumes specific values. From 
Eq. ( |Tl"| ) it follows that fj,L 2 /(2n) 2 6 N is a necessary 
condition for a disconituity in the particle number N. 
For ixL 2 /(2tt) 2 <£ N, N is constant since the number of 
momentum modes enclosed by the Fermi sphere remains 
constant. However, we observe that there exists a subset 
of integer values for which no jump in the particle num- 
ber N occurs, e. g. ^L 2 /(2n) 2 = 7 and fiL 2 /(2n) 2 = 15. 
For these values of nL 2 / (2ir) 2 , no 3-tuple (ni,n2,ris) ex- 
ists such that the sum n 2 = n\ + n\ + n\ is identical 
to these values. Thus, the particle number remains con- 
stant. In any case, we refer to a jump in the particle 
number (density) as a shell effect in the following. 

Let us now discuss the dependence of universal quan- 
tities O/cf on J and L, where O has the dimension of 
energy. In this section, we only discuss finite-size and 
particle-number effects on the Bertsch parameter. The 
corresponding effects on the fermion gap will be discussed 
in the next section. We shall see that the gap shows a 
stronger dependence on the volume size, depending on 
the actual value of J. In Fig. [4] we present our results 
for the (normalized) Bertsch parameter £/£oo as a func- 
tion of the (average) particle number N and as a func- 
tion [il? /(2tt) 2 for various values of the dimensionless 
source ji, = JL 7 / 2 . Here, ^ denotes the Bertsch pa- 
rameter for J — > in the thermodynamic limit. In order 



For the fit, we have used 25 equidistant values of Snj be- 
tween Snj Ri 0.006 and Snj fti 0.015. 

Here, we have used 10 equidistant values of Snj between Snj RJ 
0.006 and Snj 0.011 for the fit. 



For a free Fermi gas at vanishing temperature, N only assumes 
integer values. 
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FIG. 4: Top panel: Normalized Bertsch parame- 
ter £(jL,N)/£oo as a function of the particle number N 
for j L = JL 7 /2 = 0.1,0.25,0.5,0.75,1.0 x 10 7/2 (from top 
to bottom) for a fixed spatial extent L of the cubic vol- 
ume, L = Lo = 10(c/eV). Note that different values of N 
correspond to theories with different densities. Bottom panel: 
Normalized Bertsch parameter ^(Jl, /-i£o)/£°° as a function 
of /zLo7(27r) 2 for the same values of JL 7 J 2 . 



FIG. 5: Top panel: Normalized Bertsch parame- 
ter £(jL,N)/£oo as a function of the particle number N 
for L = 5, 7.5, lOc/eV (from bottom to top) for a fixed value 
of the dimensionless source jl = JL r/2 = 0.1 x 10 7/2 w 316. 
Bottom panel: Normalized Bertsch parameter £(jL,N)/£oc 
as a function of nL 2 /(2n) 2 for L = 5, 7.5, lOc/eV (from bot- 
tom to top) for a fixed value of the dimensionless source ji, = 
JL 7/2 « 316. 



to obtain these results we have fixed the spatial extent L 
of the volume, L = Lq. In terms of absolute values, 
we have chosen L = lOc/eV ~ 1.97 x 10 -7 m ^> ob, 
where as ~ 5.3 x 10~ n m is the Bohr radius, see Ref. 
for our conventions. This choice for Lq fixes the scale in 
our study and implies that we measure \i as well as J in 
units determined by the length scale L . The density is 
given by n = N/Lq. 

Since we keep J and L fixed in Fig. [4j our results for 
the Bertsch parameter depend solely on fi. Thus, an 
increase of the particle number N corresponds to an in- 
crease in the chemical potential fi. For large values of (i, 
we have J/i~ 7 / 4 ~ (/xL 2 ,) -7 / 4 — > for a fixed source J 
and L = L Q . Thus, we find AT) — > ^ for \iL\ > 1. 
For small values of N, shell effects (i. e. discontinuities) 
are clearly visible in our results for the Bertsch param- 
eter. For increasing particle number N (and fixed Lo), 
however, we find that shell effects are washed out and the 
Bertsch parameter converges rapidly to its value in the 
thermodynamic limit. For the studied values of JL Q ' 
and N > 50 (jjtl% / \2tt) 2 > 2), the relative shift 5£ of the 
Bertsch parameter behaves as 12 



(16) 



Note that k F L ~ TV 1 / 3 . 



where cn < and c M < are constants. The large-A^ 
behavior of 5£ follows immediately from the definition of 
the Bertsch parameter, see Eq. (12). Moreover, we find 
that £ (jL , N) /^oo is already close to its value in the ther- 
modynamic limit (£0l, A0/£x> > 0.98), if fiL 2 /(2Tr) > 5 
(k F L > 18), i. e. N > 200 for the studied values of the 
source J. In the next section we shall discuss how order- 
parameter fluctuations alter these predictions. In any 
case, the convergence behavior depends on J. However, 
this dependence appears to be weak for the values of J 
considered in this work. 

For a fixed value of N < 200, we observe a clearly 
visible dependence of £(.7l,-W)/£oo on the source J. In 
fact, £(Jl> N)/£oo becomes smaller for increasing J. This 
behavior is compatible with our observations in the ther- 
modynamic limit where we have also found that £ be- 
comes smaller for increasing values of the source term. 
For finite values of N, we find that the dependence of the 
Bertsch parameter with J depends on the actual value 
of [xL\. In fact, this dependence becomes stronger for 
smaller values of [xL\. Note that a fixed value of fiL\ 
does not correspond to a fixed (average) particle num- 
ber N. This can be readily seen from Fig. |4} In the 
bottom panel, the actual positions of the discontinuities 
in £(j'l, M-koVCoo do not depend on J. On the other 
hand, the positions of the discontinuites in M-^o)/^<» 
depend on the particle number N (L e. on the den- 
sity n — N/Lq), see top panel of Fig. E] 

In Fig. [5] we show £/£oo as a function of N as well 
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FIG. 6: (Normalized) Bertsch parameter as a function 
of k^L for a fixed density n. The blue symbols depict the 
results for j = Jk~ 7/2 w 0.009 with fc P = 2eV/c (n w 
35.2 x 10 12 /cm 3 ), whereas the black symbols depict the results 
for j = Jk~ 7/2 = 0.5 with k F = leV/c (n w 4.4 x 10 12 /cm 3 ). 
Open and filled symbols correspond to the mean-field approx- 
imation and to our approximation including order-parameter 
fluctuations, respectively. 



as of <iL 2 for L = L /2, 3L /4, L , where the value of 
the dimensionless source j = JL 1 1 2 has been kept fixed. 
We find again that the results become independent of 
the actual value of j for large N. For a small fixed 
value of fj,L 2 (~ small iV), we observe that the Bertsch 
parameter decreases when the volume size is decreased. 
This is similar to the behavior of £ found in the infinite- 
volume limit, see Eq. ( fl4| ). Keeping \iL 2 as well as JL 1 / 2 
fixed, we find that the (average) particle number in- 
creases when we decrease the size of the volume. Since 
we have £ = \iL 2 /(3n 2 N) 2 / 3 , it follows that the Bertsch 
parameter becomes smaller in smaller volumes, provided 
that we keep \xL 2 fixed. 

Finally, we would like to discuss the behavior of the 
Bertsch parameter as a function of the dimensionless 
quantity kpL for a fixed density. In Fig. [6] we show our 
corresponding results for ^(jjN)/^ for various values 

— 7/2 

of J, where j = Jk F ■ Filled circles depict our mean- 
field results. The results from an approximation beyond 
the mean-field limit (open circles) will be discussed in the 
next section. For small values of kpL, finite-size effects 
are clearly visible. In this regime, the behavior of £/£oo 
as a function of kpL is dominated by shell effects. As dis- 
cussed above, these shell effects are washed out for large 
values of the source. For large values of k^L, the Fermi 
gas approaches the thermodynamic limit, as it should be 
for a fixed density. This limit corresponds to a large- N 
limit since k^L ~ TV 1 / 3 . Note the difference to Fig. |ZJ 
where we have fixed the size of the volume rather than the 
density. In the present case, the Bertsch parameter does 
not approach its value for J — > in the thermodynamic 
limit for large values of k^L. This is due to the fact that £ 



FIG. 7: (Normalized) Dimensionless fermion gap as a func- 
tion of kpL for a fixed density n. The blue symbols de- 
pict the results for Jk^ 7 ^ 2 ~ 0.009 with fcp = 2eV/c (n « 
35.2 x 10 12 /cm 3 ), whereas the black symbols depict the results 
for j = Jk~ 7/2 = 0.5 with fc P = leV/c (n w 4.4 x 10 12 /cm 3 ). 
Open and filled symbols correspond to the mean-field approx- 
imation and to our approximation including order-parameter 
fluctuations, respectively. We add that, in the mean-field ap- 
proximation, Jfcp 7 / 2 w 0.009 corresponds to Snj w 0.006 in 
the infinite- volume limit, whereas j = Jk^ 7 ^ 2 = 0.5 corre- 
sponds to 8nj w 0.154 in the infinite- volume limit. These 
different values for Snj then translate into different values for 
the fermion gap in the thermodynamic limit, see Fig. [2] 

depends on J for a fixed density, even in the limit L — > oo, 
see Fig. [I] As in our study of £{jl,N) /£oo as a function 
of N, we observe that the value of k-pL above which £ has 
effectively assumed its value in the thermodynamic limit 
depends slightly on the value of the source J. In any case, 
we conclude that finite-size and particle number effects 
on the Bertsch parameter are rather weak. 



IV. BEYOND MEAN-FIELD AND THE ROLE 
OF FLUCTUATIONS IN A FINITE VOLUME 

In this section, we discuss the effects of corrections be- 
yond the mean-field approximation. These corrections 
are associated with fluctuations of the Nambu-Goldstone 
bosons. The latter arise due to the spontaneous break- 
down of the U(l) symmetry and dominate the IR physics. 
We shall see that such fluctuations play a prominent 
role in a study of finite-size and particle-number effects. 
Whereas (4>}j—o is not necessarily zero in a finite vol- 
ume in the mean-field approximation due to the absence 
of Nambu-Goldstone fluctuations, we expect (<f>)j — » 
for J — > for any finite value of L, see our discussion 
in the previous section. This behavior can be indeed 
observed in our RG analysis once order-parameter fluc- 
tuations are taken into account, cf. also previous RG 
studies of finite-volume effects in QCD low-energy mod- 
els g3HH]. 
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Before we discuss our numerical results in detail, let 
us briefly comment on the approximation which under- 
lies our studies in this section. Here, we do not discuss 
the derivation of the flow equations. Details concerning 
the derivation of the flow equations as well as explicit ex- 
pressions thereof can be found in Refs. [24l l25l l27 l l29l I3T] . 
In this work, we have simply generalized the flow equa- 
tions given in Ref. [29 to the case of a finite source J 
and a finite volume V = L 3 . To be more specific, we 
include the bosonic loops in the RG flow of the order- 
parameter potential. Moreover, we take into account 
the running of the wave-function renormalization see 
our discussion above. Such an approximation is sufficient 
for an initial RG study of the effect of order-parameter 
fluctuations on physical observables in a finite volume. 
Analogously to our mean-field analysis, finite- volume ef- 
fects are taken into account by replacing the (continu- 
ous) squared spatial (loop) momenta q 2 with discrete 
momenta q 2 = {2it) 2 n 2 / L 2 . The source J is again in- 
cluded in the RG flow on all scales k. This renders the 
mass mo of Nambu-Goldstone bosons finite and yields 
an additional contribution to the mass mjt of the radial 
mode: 

m G = 4?f= i m R = m G + 2 KP0 ■ (17) 

V z Po 

As in our mean-field study, the fermion gap has only an 
implicit dependence on the source J. The explicit form of 
the flow equations for the order-parameter potential and 
the bosonic wave-function renormalization can be found 
in App. [£] 

We begin our discussion of the effect of dynamical 
Nambu-Goldstone bosons with an analysis of the (nor- 
malized) Bertsch parameter £(j, kpL)/^ as a function 

—7/2 . rn 

of kpL, where j = Jk F , see Fig. 6 In the thermo- 
dynamic limit and for J — > 0, we rind « 0.55 in 
the present approximation. We observe that our results 
for £(j, fcF-k)/£oo are in a good agreement with the cor- 
responding results from the mean-field approximation. 
For the presently studied values of J and k-p we find that 
finite- volume effects are clearly visible for kpL < 8. For a 
given value of J and kpL > 8 (N > 20) the (normalized) 
Bertsch parameter ^{j^kpL)/^ then converges rapidly 
to its value in the thermodynamic limit. 

While our present results for the relative shift of the 
Bertsch parameter in a finite volume are in reasonable 
agreement with our mean-field results, we find that the 
fermion gap is more sensitive to the inclusion of order- 
parameter fluctuations. In Fig. [7] we show our results 
for the fermion gap as a function of kpL for a fixed den- 
sity. Filled circles depict the mean-field results. The 
results from our study beyond the mean-field limit are 
depicted by open circles. First of all, we observe shell 
effects in the fermion gap as well. However, these effects 
are less pronounced than in the case of the Bertsch pa- 
rameter. This can be understood from the fact that the 
Bertsch paramter is directly related to the density, while 
the fermion gap depends only implicitly on the density 



of the gap. In the limit kpL ^> 1 our results for the 
fermion gap converge to their values in the thermody- 
namic limit, as it should be. However, the fermion gap is 
different for different values of J, see Fig. [2j For J — » 
and kpL — > oo, we find Aoo = A/ep ~ 0.58 in the 
present approximation. For large values of the dimen- 
sionless source Jk F , we observe that the gap is almost 
independent of kpL. In this regime, the wave-length as- 
sociated with the lightest excitations in the spectrum, 
namely the Nambu-Goldstone bosons, is (much) smaller 
than the extent L of the volume, i. e. l/mg < L. Thus, 
fluctuations of the Nambu-Goldstone particles are not 
strongly affected by the presence of a boundary. In fact, 
we observe that corrections beyond the mean-field ap- 
proximation only affect the absolute value of the fermion 
gap in this regime, whereas the qualitative dependence 
on kpL remains unchanged compared to the mean-field 

7/2 

limit. Lowering Jk F , we find that the fermion gap 
already starts to deviate from its value in the thermody- 
namic limit for large values of kpL. To be specific, we find 
that finite- volume effects are clearly visible for kpL < 15 
(N < 110), see black symbols in Fig. [7] We add that 
shrinking the volume size roughly corresponds to an in- 
crease of the temperature in the Euclidean formalism. 
Therefore the decrease of the gap for small values of kpL 
can be viewed as a "melting" of the condensate (</>) j. 

Finally, we would like to discuss the approach to 
the thermodynamic limit for a situation in which the 
limit kpL might not be easily accessible numerically, e. g. 
in lattice simulations. In principle, the results in Fig. [7] 
suggest two ways to approach the thermodynamic limit. 
Apart from the direct approach of simplifying increas- 
ing kpL for a fixed density, a second approach is opened 
up by studying the Fermi gas with large values of the 
source J, such that 1/tog "C L. The results for the 
fermion gap are then essentially identical to their values 
in the limit kpL — ¥ oo. From a (linear) fit to the (dimen- 
sionless) fermion gaps obtained in this way, we can then 
extract the (dimcnsionless) fermion gap for J — > 0, see 
our discussion of Fig. [2j 

From our present results together with the observed 
insensitivity of the mean field approximation on the size 
of the system, we see that the finite-size effects are dom- 
inated by the bosonic order parameter fluctuations. To 
conclude this section, we therefore comment on the re- 
liability of the truncation of the bosonic sector. Here, 
we have only included leading-order corrections in the 
derivative expansion. However, we expect that correc- 
tions from higher orders in the derivative expansion only 
affect our results for the relative shift of the Bertsch pa- 
rameter and of the fermion gap on a moderate quanti- 
tative level but not qualitatively, since standard power 
counting arguments should be reliable for this system 
away from criticality. However, we would like to stress 
that our RG approach is by no means restricted to a 
low-order derivative expansion of the effective action. 
In the context of ultracold gases, the full momentum- 
dependence of propagators and vertices has been resolved 
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in, e. g., Refs. [571 [731 [75] , where a vertex expansion of the 
effective action has been studied. Apart from dropping 
higher orders in the derivative expansion, we have not 
considered the possibility of an inhomogeneous ground 
state in our study. This would indeed require the inclu- 
sion of higher orders in the derivative expansion. We ex- 
pect that corrections from inhomogeneities are sublead- 
ing, at least for reasonably large values of N. For small 
particle numbers N, on the other hand, inhomogeneities 
in the ground state may play an important role. Despite 
these shortcomings of our analysis, we believe that our 
present study allows us to explore the onset of finite-size 
and particle number effects in Fermi gases and therefore 
our results may help to understand better finite-size ef- 
fects in lattice (Monte-Carlo) simulations. 



V. CONCLUSIONS AND OUTLOOK 

In this work we have analyzed finite-size and particle- 
number effects in a Fermi gas confined in a periodic 
box V — L 3 . We have computed the dependence of the 
Bertsch parameter and the fermion gap on the average 
particle number N and the volume size V. Moreover, 
we have studied the dependence of these observables on 
an external pairing source J in a finite volume as well 
as in the infinite-volume limit. In the present analysis, 
we have introduced such a source to control symmetry 
breaking in a finite volume. For our studies, we have em- 
ployed non-perturbative RG techniques which allow us to 
go systematically beyond the mean-field approximation 
by means of a derivative expansion of the effective action. 
This opens up the possibility to study the effects of order- 
parameter fluctuations on physical observables, such as 
the Bertsch parameter and the fermion gap, in a system- 
atic fashion. We find that there are only mild corrections 
to the volume dependence of the Bertsch parameter aris- 
ing from the inclusion of such effects, whereas the behav- 
ior of the fermion gap as a function of kpL is strongly 
affected by order parameter fluctuations. To be specific, 
the Bertsch parameter is already close to its value in 
the continuum for kpL > 8 (N > 20) for a given fixed 
density. The observed rapid convergence of the Bertsch 
parameter for N > 20 is in accordance with recent results 
from Quantum Monte-Carlo simulations of a resonantly 
interacting Fermi gas in a periodic box [33] . On the other 
hand, we find a significant dependence of the fermion gap 
on k-pL for a fixed density, which becomes stronger for 
small values of the external source. The observed "melt- 
ing" of the fermion gap for decreasing k-pL is associated 
with strong fluctuations of the Nambu-Goldstone bosons. 

Alternatively to a study of the fermion gap Aj ~ (4>)j, 
which plays the role of an order parameter, one might be 
interested in a study of in the absence of a source J. 
The latter expectation value is not necessarily zero in a fi- 
nite volume for J = 0. Nonetheless, finite-volume effects 
are also visible in (\<f>\). In fact, it has been found in, 
e. g., lattice studies of relativistic 0(2) models that (\<j>\) 



vanishes for small volume sizes in the same way as it 
vanishes for high temperatures, see e. g. Ref. [76]. De- 
pending on the setup of a given lattice simulation of an 
ultracold Fermi g detailed analysis of the onset of 
finite- volume effects in (\<f>\) might therefore be worth- 
while as well. A corresponding RG study is deferred to 
future work. 

Finally, we add that the observed volume dependence 
of the fermion gap may affect the determination of the 
critical temperature in the unitary limit in Monte-Carlo 
simulations. The present work aims to set the stage of 
a detailed RG scaling analysis of the critical behavior of 
a resonantly interacting Fermi gas, including a study of 
the dependence of the critical temperature on the vol- 
ume size and the particle number. Another interesting 
extension of the present RG approach represents a study 
of a Fermi gas based on the canonical ensemble rather 
than the grand canonical ensemble which we have con- 
sidered here. This might be possible in a similar way as 
done in lattice QCD simulations, see e. g. Ref. [77]. Such 
an analysis may then provide additional insights into the 
many-body physics of ultracold atomic gases and may 
help to further bridge the gap between Monte-Carlo sim- 
ulations and continuum approaches. 
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Appendix A: RG Equations Beyond the Mean-Field 

Limit 

Beyond the mean-field approximation, the RG flow of 
the order-parameter potential U receives contributions 
from a purely bosonic loop. To be specific, the flow equa- 
tion reads 



d t U{p,J,L,ii) = n vP U' -2k 5 (B> +B<)s F 

~5 



+ k 5 (l- 2 ^)B B s B , (Al) 



where p = Ztpp, rj v = —dt In z\ and 



(kLy 



3tt 2 



This function counts the (bosonic) momentum modes 
and corresponds to the mode-counting functions for the 
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fermions defined in Eqs. |8]) and ([£]). Recall that we use 
periodic boundary conditions for the fields. The func- 
tion Sb is given by 



see Ref. The RG running of Zjj, is determined by 

the anomalous dimension r\ v : 



•SB 



1 + 0J1 
1 + UJ 2 



1 + C02 
l + U)i 



(A2) 



where wi = U'/k 2 and uj 2 = (V + 2pU")/k 2 ; the prime 
denotes the derivative with respect to p. The quan- 
tities uj\ and W2 are related to the (bare) masses mc 
and ?7iR which we have defined in Eq. (|17|: 



127T 2 (1+W 3 )t 



PS 



B<)(2-u 3 ), (A5) 



Po 



m G 

z\k 2 



Here, we have used that 

,2 / 



w 2 = — 



17) 



(A3) 



where W3 = / 
defined in Eqs 



^po/fc 4 - The functions Bp and ±Sp are 
{[sj) and (|9]), respectively. Note that a non- 
trivial RG running of Zip induces a non-trivial RG-scale 
dependence of the (dimensionless) Yukawa coupling h 2 = 



(A4) h 2 /(Z!pk), see our discussion in Sect. 
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